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Abstract 

We investigate the classical scattering dynamics of the driven elliptical billiard. Two funda- 
mental scattering mechanisms are identified and employed to understand the rich behavior of the 
escape rate. A long-time algebraic decay which can be tuned by varying the driving amplitude is 
established. Pulsed escape rates and decelerated escaping particles are generic properties of the 
harmonically breathing billiard. This suggests time-dependent billiards as prototype systems to 
study the nonequilibrium evolution of classical ensembles encountering a multitude of scattering 
processes off driven targets. 

PACS numbers: 05.45.-a,05.45.Ac,05.45.Pq 
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Many models of statistical mechanics can be reduced to billiard-type dynamical systems. 
Seminal results like the justification of a probabilistic approach to statistical mechanics {l|, 2| 
or a connection between the escape from a billiard and the famous Riemann hypothesis 3| 
were obtained. Modern studies of billiards are, among others, motivated by corresponding 



experiments including ultracold atoms confined in a laser potential 



billiards 



mm 



microwave 



or mesoscopic quantum dots Equally for the design of directional 



micro-lasers, billiards are of immediate relevance 



12l |. From the theoretical point of view. 



investigations on the classical and quantum properties of billiards have pioneered the fields 
of quantum chaos, modern semiclassics and transport at the mesoscopic scale (see js] and 
refs. therein). 

The escape rate, being the fraction of remaining particles as a function of time, is a 
characteristic of open billiards which is both experimentally accessible and important for 
transport properties This key property allows to probe the dynamics from the outside 



151] • Integrable 



and has been studied thoroughly for billiards with a static boundary 
systems exhibit an algebraic, while fully chaotic billiards show an exponential decay of the 
escape rate 13|] (altough it is known that in special cases, like the Bunimovich stadium, 
there is an algebraic long-time behavior in the decay, due to the slow transport of particles 
close to the marginally stable bouncing ball orbits, cf. e.g. Q). Here we focus on the case 
where the billiard shape changes in time according to a certain law. Such a driven billiard 
not only leads to a higher dimensional phase space of the scattering processes but also to 
a non-conservation i.e. time-evolution of the energy. Driven billiards represent prototype 
systems for the evolution of ensembles of particles in a closed driven environment where 
multiple scattering off the driven boundary takes place thereby leading to a dynamical non- 
equilibrium state. Extremely little is known on the properties of suc h sy stems, the few 



existing investigations dealing with aspects such as Fermi acceleration [16 
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19| and 



22l |. In spite of the above- 



principal structures of the corresponding phase space [20|, 
motivated interest escape rates of driven billiards have not been addressed up to date. Using 
atom-optical techniques such as acousto-optic deflectors or even in the case of microwave 
cavities moving boundaries are well within reach of experiment. 

We investigate the classical scattering dynamics and the time-evolution of ensembles of 
particles in a harmonically driven elliptical billiard. The decay of the escape rate is traced 
back to the underlying scattering mechanism by identifying two fundamental scattering pro- 
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cesses being key ingredients for the time-evolution. The escape rate behaves asymptotically 
as Nc{t) ~ t'"^^ and it is shown that the decay constant wc can be changed continuously 
by varying the driving amplitude C. Moreover, phenomena such as a pulsed escape, demon- 
strate the richness of the properties of driven billiards and suggest specifically the driven 
ellipse as a prototype system for the nonequilibrium evolution of ensembles that experience 
multiple scattering processes with moving targets. 

As a precursor for the driven system let us discuss some relevant features of 
the static elliptical billiard. The boundary B of an ellipse is given hj B = 
{{x{ip) = Acosip,y{ip) = B sirup) \0 < (f < 27i} with A > B > 0, thus A and B being the 
long and the short half-diameter, respectively. The dynamics in the ellipse is completely 
integrable 20|], see the Poincare surface of section (PSS) shown in Fig. [H In addition to 
the energy, there is a second constant of motion 



p^(l + (l-e^)cotV)-e^ 
"^^^'^^^ l + (l-.^)cotV-^^ ' 



restricting the orbits to invariant curves in phase space, where e = a/1 — B"^ / A^ is the 
eccentricity and p = cos a {a is the angle between the tangent at the boundary and the 
trajectory of the particle). Ffcp, p) can be interpreted as the product of the angular momenta 
(PAM) about the two foci [23]. There are two different types of orbits in the ellipse, librators 
and rotators, see Fig. [H Librators cross the x-axis between the two foci and touch repeatedly 
a confocal hyperbola. They possess values of F between — 6^/(1 — e^) and zero and their 
motion is restricted to a limited range of (f. Rotators travel around the ellipse, eventually 
exploring every value of ip, repeatedly touching a confocal ellipse. They have values of F 
between zero and one. 



ESCAPE RATE 

Let us consider the escape rate of the static billiard by placing a small hole at the very 
right {(p = 0) of the ellipse {A = 2, B = lin arbritrary units) and iterating the corresponding 
2D discrete mapping (for the phase space variables ip and p) numerically. We employ an 
ensemble of 10^ particles with initial conditions being uniform randomly distributed in ip, a- 
space. The result is shown in Fig. [2] (curve with C = 0.00). The decay approaches a 
saturation value Ns{e) which is caused by particles traveling on librator orbits that are 
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not connected with the hole j^. For e = the saturation value is zero and it increases 
monotonically with increasing e. The analytical dependence of Ns on e can be obtained by 
determining the number of librator type particles starting from the region bounded by the 



separatrix 



NJe] 



2tt 



dip arccos 



1 + (1 - £2) cot^ 



(2) 



This expression is in excellent agreement with the results of the numerical simulations. One 
can therefore conclude that varying e allows us to control the total number of particles being 
emitted. 




FIG. 1: Phase space of the ellipse, the invariant curves are the contour lines of F(ip,p) (upper 
part) and typical trajectories in coordinate space (lower part). 



We now apply harmonic oscillations to the boundary of the ellipse 



{Aq + C sin(ct;t + 5)) cos (f 
{Bq + C sin(co't + 5)) sin 



(3) 



C > is the driving amplitude and 5 is a phase shift which will be set to zero in the following. 



The dynamics of the particles is now governed by a 4D discrete mapping 20|, l2l(| : the time 
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and position of a particle-boundary collision is described by the pair (yj, t, the energy and 
direction of motion hj v = {vx,Vy). The time tn+i for the n + 1th collision is determined 
implicitly by 



+ C sm{ut + 5) J \Bq + C siTi{ut + 5) ) ' 

where the smallest tn+i > tn that solves (jlj) has to be taken and Xn = {xn,yn) is the nth 
collision point. The n+1 collision point is given by Xn+i = Xn + Vn(tn+i — tn) and ipn+i can 
be obtained by inverting ([3]). The velocity immediately after the collision is described 
by Vn+i = f„ — 2 [hn+i ■ {vn — "U-n+i)] " fi-n.+i, wherc Un+i IS the boundary velocity and n„+i 
the inward pointing normal vector of the coUisional event occurring at time tn+i and position 
Xn+i- When iterating the underlying implicit mapping numerically the determination of 
from eq. (HI) requires sophisticated numerical techniques due to many neighboring roots. 
As a result the corresponding simulations are computationally very demanding. 



We focus on the escape rates of monoenergetic ensembles consisting of 10^ particles ^ 
with uj = 1, Aq = 2 and Bq = 1 for different values of the driving amplitude C. Two 
relevant cases can be distinguished: |i;o| ~ ^^C being the intermediate velocity ensemble 
(IVE) for which the initial particle velocity is of the order of the maximum velocity of the 
billiard boundary and \vo\ ^ uC being the high velocity ensemble (HVE). Naturally, it 
would be also interesting to examine the case |i^o| ^ ^C. However, the first few collisions 
then accelerate the particles to velocities ^ uC and, after a short time, we encounter the 
situation of the IVE. Nc{t) of the IVE and HVE are shown in Figs. [2] and [3l respectively. 
In both cases, there exists no saturation value as observed in the static case. A short-time 
exponential decay is followed by a transient and the long-time behavior {t > 10^) is to a 
good approximation algebraic, Nc{t) ~ t^"^^' , at least for the case of the IVE (we remark 
that this algebraic decay has been numerically shown to exist for much longer times than 
illustrated in Fig. [2]). The decay constant wc increases monotonically with increasing C 
(this fact is based not only on the four values of the drivig amplitude C shown here, but 
on simulations carried out for 20 values of C between 0.01 and 0.30). It can be therefore 
concluded that tuning the driving amplitude allows us to control the decay constant wc in 
the long-time behavior of the decay. The division of the behavior of the escape rate into 
the above-mentioned regimes is even more pronounced in the case of the HVE, see Fig. [3l 
although the long time tail shows substantial deviations from an algebraic behavior. 
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FIG. 2: Fraction of remaining particles Nc{t) of the IVE as a function of time for different values 
of the driving amplitude C. The inset provides a semi-logarithmic plot. 

The saturation value Ns{e) observed in the static case is caused by librating orbits which 
are gradually destroyed by the driving ([3]) thereby resulting in a non- vanishing rate Nc{t) 
even for large times. Rather than trying to examine the 4D phase space (since 2D inter- 
sections or projections of this 4D space are nonrepresentative and difficult to interpret), let 
us analyze the destruction of the librator orbits as follows. Orbits (v^i,^^) in the driven 
ellipse can be compared to the corresponding ones of the static ellipse {'^i,Pi) by inspect- 
ing the angular momentum F{(p,p). In contrast to the case of the static ellipse where 
F{ipi,pi) = const. Wi, F{ip[,p'^ ^ F{ip'-,p'j) {i ^ j) in the driven case, i.e. F is no longer 
a constant of motion. The difference AF = Fj+i — Fi (see Fig. [1]) upon a collision is a 
measure of whether a librator orbit approaches {AF > 0) the separatrix or whether it moves 
in phase space towards (AF < 0) the corresponding elliptic fixed points of the static system. 
Two fundamental scattering processes can be identified that destroy librators by changing 
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FIG. 3: Fraction of remaining particles Nc{t) for the HVE, 27r-oscillations of the decay are shown 
in the inset. 

them into rotator orbits: 

1. Vertical process: The angle of incidence of a particle upon collision does not coincide 
with the reflection angle due to momentum momentum transfer by the moving billiard 
boundary. This momentum change Ap causes a vertical displacement of the particle 
in the PSS. 

2. Horizontal process: A particle that would hit the boundary at a certain point (f in 
case of a static billiard, hits the boundary in the driven case at (/?', simply because 
it has moved. Here no change of the momentum takes place. This corresponds to a 
horizontal displacement A(f of the particle in the PSS. 

These processes are fundamental in the sense that every change AF can be decomposed into 
AF — AFh + AF„, where AFh^^ denote the individual changes caused by the horizontal, 



the vertical process, respectively. 

In general, these effects do not appear isolated but a combination (Ay?, Ap) of both 
will occur at a single collision. The value AF/j ^, at a single collision can be calculated 
exactly [24]. The changes AF^^v increase monotonically with the driving amplitude C. 
Additionally, AF^ depends monotonically on the driving frequency u. It is important to 
note that collisions occurring during the expansion/contraction period of the ellipse always 
lead to AFy > 0/ AF^ < 0, whereas such a clear distinction is not possible for AFh. 

For a single particle undergoing many scattering processes, the effective change AFe 
after a certain time depends on the sequence of these processes, hence AFe = AFi + AF2 + 
• ■ ■ + AFn after n collisions. When regarding an ensemble of particles, the effective changes 
( AFe)j (where the index j indicates the jth particle) after n collisions can vary significantly 
from particle to particle, since the sequence of the individual AFi will vary substantially for 
different initial conditions which is due to the nonlinear dynamics of the underlying discrete 
mapping. 

With the just presented considerations, we can explain qualitatively the observed escape 
rates and especially the disappearance of the saturation value. Let us focus first on the HVE. 
The initial fast decay {t < 5) is due to the rotator orbits that are connected with the hole 
(in the PSS) and escape very rapidly. The long-time decay (t > 10) is caused by particles 
starting on librator orbits that have been scattered through horizontal and vertical processes 
across the separatrix. The closer a particles initial orbit lies near the elliptic fixed points (of 
the static billiard), the longer it takes until the effective change AF is large enough to reach 
the separatrix {F = 0), see Fig. HI The monotonic dependence of the AFi on C explains 
the increasing emission rate Nc{t) with increasing C, since at a given time t, the number 
of particles that can participate in the decay is larger for larger values of C. The decay in 
the transient region (5 < t < 10) is caused by a superposition of the tail of the initial fast 
exponential decay and the onset of the slow algebraic decay. With the same arguments, the 
decay of the IVE can be explained qualitatively. Since O{vo) ^ 0{max{u)), the resulting 
changes AF are much larger for the case of the IVE compared to the ones of the HVE. This 
leads to a very early onset of the slow algebraic decay and consequently the transient region 
is broadened. 

The above presented considerations are confirmed by the results shown in the right inset 
of Fig. m where exemplarily the average escape time and the standard deviation a of the 
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IVE (C = 0.10) as a function of the initial PAM F{(po,po) are plotted. Particles with values 
of F{(Pq,pq) between zero and one (rotators) have approximately the same average escape 
times < tesc >, but the large (compared to < tesc >) standard deviation a indicates that 
the individual escape times can vary significantly from particle to particle. The average 
escape time of particles starting on librator orbits {F{(fQ,pQ) < 0) increases with decreasing 
F{(po,Po), since the required AFe to change a librator into a rotator becomes larger and 
more and more collisions are necessary to reach these values of AFe. 

In the inset of Fig. [3l a modulation of the escape rate with period T = 2tx can be 
seen, being exactly the period of the applied driving law (j3]). Specifically, for t > 10, 
where all particles starting on rotator orbits have already escaped, Nc{t) ~ const, during 
approximately 11/12 (empirically observed) of one period and subsequently Nc{t) ^ 
during a time interval T/12 only. From this behavior it is evident that the ellipse operates 
from a certain time on as a pulsed source of particles. These repeated intervals are centered 
around points tm of maximal extension of the ellipse, tm = (4m + l)vr/2, m = 2,3,4, .. . 
During the expansion period, dominantly vertical but also horizontal processes turn librators 
into rotators. The moving ellipse remains for a comparatively long time period in the vicinity 
of the extremal configuration at and consequently the newly created rotators escape. 
Therefore, the dynamics is effectively probed during these short time intervals centered 
around tm- During the contraction period, the librators are stabilized via vertical processes, 
consequently Nc{t) ~ during 11/12 of a period T. 

An astonishing feature of the escape of the HVE is displayed in Fig. HI where the distri- 
bution p(|f I) of the escape velocities and the correlation of the escape time vs. the escape 
velocity (left inset) after a certain time is shown for C = 0.25. There is a strong asymme- 
try of p(li'l) around the initial velocity \vq\ = 100. Particularly, particles with tesc > 10, 
being originally exclusively librators, possess escape velocities \v\ < \vq\. The origin of this 
behavior are vertical processes which turn librators into rotators that subsequently escape. 
The latter process however takes place during the expandin g ph ase of the moving ellipse 
which leads exclusively to a decrease i.e. loss of energy 25|, |26|]. The escape velocities of 
particles starting on rotator orbits lie on a serpentine band, whose oscillations match the 
periodic oscillations of the ellipse. We therefore conclude that the breathing ellipse acts as 
a decelerator for particles in the HVE regime. 

In conclusion, we have explored the dynamics of a driven elliptical billiard with a focus 
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FIG. 4: Distribution p{\v\) of the velocity of the escaped ensemble (HVE, C = 0.25). Left inset: 
the escape time vs. the escape velocity is plotted. Right inset: average escape time of the IVE 
(C = 0.10) as a function of the initial F{ipo,po). 

on the scattering mechanisms and their impact on the escape rate. Observed phenomena 
such as the long-time algebraic decay which can be tuned by varying the driving amplitude, 
a pulsed escape rate and lowered velocities of the escaping particles are all caused by projec- 
tiles emanating from librator orbits. They can be understood by means of two fundamental 
scattering processes that turn librating into rotating orbits. Consequently, these effects can 
be further enhanced and manipulated by preparing suitable initial ensembles. For example, 
a high velocity ensemble initially consisting exclusively of librators restrains particles with 
an escape time t < 10 which do not show a pulsed decay and whose escape velocities are 
symmetrically distributed around the initial velocity. It is foreseeable that the application 
of different driving laws and driving modes (e.g. area-preserving oscillations) as well as the 
preparation of suitable initial (thermal) ensembles will further advance equally our under- 
standing as well as the phenomenology of driven nonequilibrium systems. Driven billiards 
serve as prototype systems in this respect. 
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